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The measurement of the pion form factor and, more generally, of the cross section for electron-positron annihi- 
lation into hadrons through the radiative return has become an important task for high luminosity colliders such 
as the $- or B-meson factories. This quantity is crucial for predictions of the hadronic contributions to (g — 2) M , 
the anomalous magnetic moment of the muon, and to the running of the electromagnetic coupling. But the 
radiative return opens the possibility of many other physical applications. The physics potential of this method 
at high luminosity meson factories is discussed, the last upgraded version of the event generator PHOKHARA is 
presented, and future developments are highlighted. 



1. INTRODUCTION 

Electron-positron annihilation into hadrons is 
one of the basic reactions of particle physics, 
crucial for the understanding of hadronic inter- 
actions. At high energies, around the Z res- 
onance, the measurement of the inclusive cross 
section and its interpretation within perturbative 
QCD |1I2) give rise to one of the most precise 
and theoretically founded determinations of the 
strong coupling constant a s P). Also, measure- 
ments in the intermediate energy region, between 
3 GeV and 11 GeV can be used to determine a s 
and at the same time give rise to precise mea- 
surements of charm and bottom quark masses fU . 



The low energy region is crucial for predictions 
of the hadronic contributions to a lll the anoma- 
lous magnetic moment of the muon, and to the 
running of the electromagnetic coupling from its 
value at low energy up to Mz- 

Last, but not least, the investigation of the ex- 
clusive final states at large momenta allows for 
tests of our theoretical understanding of form fac- 
tors within the framework of perturbative QCD. 
Beyond the intrinsic interest in this reaction, 
these studies may provide important clues for the 
interpretation of exclusive decays of B- and D- 
mesons, a topic of evident importance for the ex- 
traction of CKM matrix elements. 
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2. HADRONIC CONTRIBUTIONS TO 

a M AND «qed 

The main uncertainty to and oqed is driven 
by their respective hadronic contributions, which 
are estimated though dispersion integrals 
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where the spectral function 

R(s)cx |(0|J AI |had,( 7 ))| 2 , (2) 

is obtained from experimental data of the reac- 
tion e + e~ — * hadrons. The most recent experi- 
mental result for [S] shows a 1.9c discrepancy 
with respect to the SM prediction for this quan- 
tity |(il7l8j . Alternatively, one can also use cur- 
rent conservation (CVC) and isospin symmetry 
to obtain R(s) from r decays. In the latter, a 
0.7cr discrepancy is found [H], which however is 
incompatible with the e + e~ based result. Unac- 
counted isospin breaking corrections due to the 
difference of the mass and width of the neutral to 
the charged p-meson could explain this discrep- 
ancy £8], leaving the e + e~ based analysis as the 
most reliable. 

3. RADIATIVE RETURN AT MESON 
FACTORIES 

The recent advent of $- and i?-meson facto- 
ries allows us to exploit the radiative return to 
explore the hadronic cross section in the whole 
energy region from threshold up to the nominal 
energy of the collider in one homogeneous data 
sample |9lll)| . The radiative suppression factor 
0(a/n) is easily compensated at these factories 
by their enormous luminosity. 




(a) (b) 



Figure 1. Leading order contributions to the reac- 
tion e+e- -> 7r+7r~7 from ISR (a) and FSR (b). 



In principle, the reaction e + e~~ — > 7 + hadrons 
receives contributions from both initial- and final- 
state radiation (Fig. [IJ, ISR and FSR respec- 




ts (d) 



Figure 2. Typical amplitudes contributing to the 
reaction e + e~ — » tt+tt^jj. For two photons 
emitted either from the electron/positron or the 
hadronic system. Only one representative is dis- 
played. 



tively. Only the former is of interest for the ra- 
diative return. A variety of methods to disentan- 
gle FSR from the ISR contribution have been de- 
scribed in detail in |9I11I12I13I14| . The first possi- 
bility is based on the fact that FSR is dominated 
by photons collinear to ir + or 7r~, ISR by pho- 
tons collinear to the beam direction. This sug- 
gests that one should consider events with pho- 
tons well separated from the charged pions and 
preferentially close to the beam. The second op- 
tion is based on the markedly different angular 
distributions of ISR and FSR, and the character- 
istic feature of their interference. Various charge- 
asymmetric distributions can be used for indepen- 
dent tests of the FSR model amplitude, with the 
forward-backward asymmetry as simplest exam- 
ple. Note however that at _B-factories the 7r + 7r~7 
final state is completely dominated by ISR. 

Higher order radiative corrections (see Figs. [21 
and[3Jl are also important for the precise extrac- 
tion of the hadronic cross section through the ra- 
diative return. Furthermore, simulation of indi- 
vidual exclusive channels: e.g. 7r + 7r~, K + K~, 
pp, 7r°7r + 7r~ , Air , 5ir, 671", mrr], KKir, KKirir, 
2K2K, KKi], in the low energy region (< 2-3 
GeV) requires a fairly detailed parametrization 
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(e) (f) 



Figure 3. Typical amplitudes describing virtual 
corrections to the reaction e + e _ — > 7r + 7r _ 7. Per- 
mutations are omitted. 



of various form factors. 

4. MONTE CARLO EVENT GENERA- 
TORS 

The proper analysis requires necessarily the 
construction of Monte Carlo event generators. 
The event generators EVA 9 and EVA47T [TS] 
were based on a leading order treatment of ISR 
and FSR, supplemented by an approximate in- 
clusion of additional collinear radiation based on 
structure functions. Subsequently, the event gen- 
erator PHOKHARA was developed |llll2H3j : 
it is based on a complete next-to-leading order 
(NLO) treatment of radiative corrections (16117) . 
In its version 2.0 it included ISR at NLO and 
FSR at LO for ir + ir~ and final states, and 

four-pion final states (without FSR) with some 
improvements with respect to the formulation de- 
scribed in (T5) . 

The most recent version of PHOKHARA, ver- 
sion 3.0 [T2], allows for the simultaneous emission 
of one photon from the initial and one photon 



from the final state. This includes in particu- 
lar the radiative return to 7r + 7r~(7) and thus the 
measurement of the (one-photon) inclusive tt + 7t~ 
cross section, an issue closely connected to the 
question of 7r + 7r~(7) contributions to a M . 

Recently, a new Monte Carlo event generator, 
EKHARA PB], has been constructed to simu- 
late the reaction e + e~ — > 7r + 7r~e + e - , a potential 
background of the radiative return in particular 
at lower energies. 

Experimental studies presented in |l9l[20l21l221 
[231 indeed demonstrate the power of the method 
and seem to indicate that a precision of one per 
cent or better is within reach. In view of this 
progress a further improvement of theoretical un- 
derstanding and Monte Carlo generators is re- 
quired. 

5. PHOTONIC CORRECTIONS TO a M 

The issue of photon radiation from the final 
states is closely connected to the question of 
photonic contributions to a^, see Fig. 0] (for re- 
lated discussions, see e.g. (24125) ). At the present 
level of precision of hadronic contributions to a^, 
roughly half to one per cent, these corrections 
start to become relevant. 




Figure 4. Photonic corrections to a M . 



Qualitatively the order of magnitude of this ef- 
fect can be estimated either by using the quark 
model with m u ~ rrid ss m s rs 180 MeV 
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(adopted to describe the lowest order contribu- 
tion 1223), or with m u « w m s w 66 MeV 
(adopted to describe the lowest order contribu- 
tion to a(Mz)) or by using tt + it~ as dominant 
intermediate hadronic state plus photons cou- 
pled according to sQED. The three estimates (see 
also j23) 

Scift (quark, 7, m q = 180 MeV) = 1.880 x 10~ 10 , 
<5a p (quark,7,m 9 = 66 MeV) = 8.577 x 10~ 10 , 
Sa^Tr+n- ,7) = 4.309 x 10" 10 , 

are comparable in magnitude and begin to be rel- 
evant at the present level of precision of ±8T0~ 10 . 
This order-of-magnitude estimate suggests that a 
more careful analysis is desirable. 

Soft photon emission is clearly described by 
the point-like pion model. Hard photon emission, 
with E 1 > 0(100 MeV), however, might be sensi- 
tive to unknown hadronic physics. Therefore the 
size of virtual, soft and hard corrections has to 
be studied separately. Let us now estimate the 
contributions from hard photon radiation to a M . 
In Fig. [3] we display the integrand 




0.5 1 1.5 2 

s(GeV 2 ) 



Figure 5. Differential contribution to ' 7 from 
7r + 7r~7 intermediate states for different cutoff 
values compared with the complete contribution 
(virtual plus real corrections, labeled 'inclusive') 
evaluated in sQED (FSR), as well as with the 
contribution from the 7r + 7r~ intermediate state. 



d CY 2 

had^cut) = K{s) i?H (S)jB cu t) (3) 

as * 6TT Z S 

for E cut = 10, 100 and 200 MeV as a function 
of y/s — m(7r + 7r~7) between the threshold and 
2 GeV. The result is compared with the com- 
plete sQED contribution, as derived from point- 
like pions, and the lowest order contribution from 
7r + 7r~. The integrated result, with E~ i between 
E cut and infinity, is displayed in Fig. |S] Contri- 
butions from the hard region, above 100 MeV, are 
clearly small with respect to the present experi- 
mental and theory-induced uncertainty. 

A measurement of 7* — * 7r + 7r~7 for variable 
\fs is desirable for an independent cross check. 
This is indeed possible with the radiative return, 
if events with two photons in the final state are 
investigated ^Sj • This was in fact one of the mo- 
tivations for extending the event generator PHO- 
KHARA to events with simultaneous ISR and 
FSR. 



6. PERSPECTIVES FOR THE RADIA- 
TIVE RETURN 

Let us now describe a number of topics for the 
development of PHOKHARA |llll2ll3| . which 
are of course intimately connected to the future 
potential of measurements through the radiative 
return. To large extent this is still a wish list, 
which will at best be partially fulfilled, depend- 
ing also on the future development in the ex- 
perimental analysis. Some of these topics are 
based on fairly straightforward extensions of the 
present framework, some would require major cal- 
culations. Inclusion of additional hadronic final 
states is only possible on the basis of parametriza- 
tions (not yet available) for the production am- 
plitudes of multi-hadronic states. The radiative 
return may, finally, give access to completely new 
phenomena. For example one may scrutinize the 
_D°_D°-system, produced in the radiative return to 
the ip" , for time-dependent mixing effects, which 
are not easily accessible in other experiments. 
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Figure 6. Integrated contribution to ajf 1 ''' as a 
function of the cutoff E cut . 



6.1. THE RADIATIVE 
MUON PAIRS 



RETURN TO 



The present version of PHOKHARA (3.0) in- 
cludes the amplitudes for 



MA* 7 i 



(4) 



in Born approximation plus those 0(a) real and 
virtual corrections which are attached to the elec- 
tron side 111112113] . Important additional correc- 
tions are expected from the second order process, 
where one hard photon is emitted from the elec- 
tron, another photon is collinear with the muon. 
These additional terms are presently being in- 
cluded in the simulation |27j . Resonances, narrow 
ones like J /if), wide ones like the p-meson and 
the vacuum polarization from virtual hadronic 
and leptonic intermediate states affects the /x-pair 
cross section. In particular the effects of the vac- 
uum polarization should and will be included in 
a detailed comparison between theory and exper- 
iment. 



6.2. BENCHMARK TESTS AND COM- 
PARISONS 

To arrive at reliable Monte Carlo programs 
suited for the experimental analysis compar- 
isons between the results of simulations, semi- 
analytical results and comparisons between pro- 
grams of different authors are of paramount im- 
portance. Inclusive semianalytical results of ISR 
NLO corrections 



are in agreement with 
PHOKHARA at a per mill level or better [T2] . 
Preliminary studies |3()| indicate agreement be- 
tween PHOKHARA 3.0 [E] and KKMC [2] for 
the p + p,~j mode at the level of two per mill 
in the region of interest, at least in the mode 
where FSR had been switched off. The agree- 
ment of partial results between PHOKHARA 
and BABAYAGA 5H| event generators in tt + tt~j 
modes is encouraging. Semi inclusive analytical 
results with specific idealized cuts have been pub- 
lished in |34j. However, no comparison with al- 
ready available numerical results of or with 
the Monte Carlo program PHOKHARA has been 
presented. Clearly, all that can only be the start- 
ing point for more detailed studies. 

6.3. HADRONIC FINAL STATES 

i) Two mesons: e.g K + K~ and K°K° 

At present only two-pion and four-pion final 
states are implemented, and the parametriza- 
tions of the form factors are optimized for Q 2 
below 3 to 4 GeV 2 . Larger values of Q 2 , up 
to 0(10 GeV 2 ) are at reach at B- meson facto- 
ries. Although the implementation of amplitudes 
and form factors for arbitrary two meson final 
states is straightforward once the generator has 
been coded for 7r + 7r~ , the remaining two-body fi- 
nal states K + K~ and K°K° are still missing, as 
well as a reliable and theoretically well founded 
formulation for the pion form factor which is ap- 
plicable both in the p and p' resonance region and 
for large Q 2 at the same time. 

ii) Three meson: e.g. 3tt, ttttt] and KKir 
Strangeness-, isospin- and CP-conservation 

alone lead to significant restrictions on the struc- 
ture of production amplitudes for three meson fi- 
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nal states. CP- and G-parity restricts the three- 
pion state (and all states with an odd number 
of pions) to the isospin zero configuration, dom- 
inated by the w-resonance and its radial excita- 
tions. (The <p as well as the J/tp contributions 
are suppressed as a consequence of Zweig's rule, 
but not strictly forbidden.) The ttttj] (and irirr]') 
mode is restricted to the isospin one channel with 
a production amplitude modeled according to the 
process 



7* — > p — > f]p (— > 7T7r) 



(5) 



Here p stands for the full set of p-meson-like ex- 
citations and the amplitude is written as 

(tttttjIJJ.O) kF p (Q 2 )F p (( P+ +p-) 2 ) x 



(G) 



The KKit mode may receive contributions from 
isospin zero and isospin one with a rich struc- 
ture of p- and K*- resonances in the KK and 
Ktt subchannels respectively. The parametriza- 
tion of amplitudes for final states with higher mul- 
tiplicities, involving K, 7T, and 7/ will quickly lead 
to enormously complicated resonance structures 
and only data-driven analysis can lead to reliable 
paramctrization of the corresponding amplitudes. 

iii) Baryons: 

The situation is simple for the measurement of 
the proton (and neutron) electric and magnetic 
form factors, 



(pp\ J^l, 0) = 7 p Fi + (p v + 

Gm = ^i + F 2 , , 
G E =F 1 + tF 2 , 



P- 



F2_ 

4m 



(7) 



with t = Q 2 j\w? B > 1. Sizeable event rates can 
be reached out to fairly high Q 2 (see Fig[7f and 
the angular distributions allow for a clean sepa- 
ration of the two form factors . 

iv) Multi-particle continuum: 

Up to now we have restricted the discussion to 
exclusive final states. The situation may become 
even more complicated for the multi-particle con- 
tinuum at larger Q 2 , say above 0(10 GeV 2 ). In 
the parton model this reaction would be modeled 
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Figure 7. Differential in Q 2 cross section of the 
process e+e~ — > ppj for two different sets of cuts. 



by e + e _ — * jqq. Hadronisation of the qq state 
could be simulated in a phenomenological way. 
However, final state radiation will be an impor- 
tant background, and elaborate methods, based 
e.g. on the markedly different angular distribu- 
tions of ISR versus FSR will have to be developed 
to analyze and separate the two contributions [see 
e -g- EE3I Eq.(3)]. In principle this method could 
also be employed to measure the R-ratio above 
charm threshold. Making use of the large boost of 
the D-mesons and their displaced decay vertices, 
the separation of charm and non-charm contribu- 
tions might be feasible, thus providing an impor- 
tant input for the determination of m c through 
sum rules From the purely technical view- 

point PHOKHARA can already now be used for 
the simulation of this reaction, by simply replac- 
ing muon mass and charge by the corresponding 
values for the charm quarks (Note, that a sizeable 
forward-backward asymmetry is predicted from 
ISR-FSR interference also in this case.). 
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6.4. NARROW RESONANCES 

The production of narrow bottonium states 
through the radiative return has been discussed in 
some detail in $5T\ (see also |38l39p . Here we shall 
concentrate on the J/^-resonance and its radial 
excitations. Neglecting for the moment radia- 
tive corrections and interference between the res- 
onance and the continuum amplitude, the cross 
section is given by 

r e Una (l + (M 2 /s) 2 /l + cos6» m \ 
° R ~ M^^\l-(M 2 /s) 2 °\\-co&B m ) 

(l- cos0 m j . (8) 

The angular region of the photon is restricted 
by | cos 6*| < cos# m and T e denotes the resonance 
partial decay rate into e + e~. For an integrated 
luminosity of lOOfb -1 this corresponds to 34- 10 4 
J/ip events, out of which 20000 events will be 
found in the channel. It is remarkable, 

that the results form the BABAR-collaboration 
for this reaction (see |21I22| ) have lead to a value 
for TeTft/Ttot which already now is more precise 
than the world average as compiled in |4U| . The 
enormous statistics gives access to a multitude of 
final states. In particular it is possible to mea- 
sure simultaneously the cross section on- and off- 
resonance. For the double ratio 

RD f = ° S (on res.)/ °* (off res.) , (9) 

one evidently expects RDy = 1, if the decay 
J/^ —* f can proceed through the virtual pho- 
ton only. This is expected for final states with 
isospin 1, like 2ir, Air, 6ir or ttttj]. Final states 
with isospin zero, like 3tt, 5tt, KK, KKtt are 
dominated by hadronic J/ip decays and, in the 
framework of perturbative QCD mediated by the 
three gluon intermediate state. 

6.5. RADIATIVE PRODUCTION 

OF ^(3770) 

A wealth of physics results has bee obtained 
from BB production at asymmetric B factories. 
The BB state is produced in a 1 — - configura- 
tion, which leads to nontrivial correlations be- 
tween the decay products of B and B mesons. 



The Lorentz boost allows to resolve the corre- 
lated time dependence of the two decays, which 
in turn is the basis for the measurement of the 
CP asymmetry. Similar studies can, in principle 
be performed for the DD system with the help of 
the radiative production of 0(3770) and its sub- 
sequent decay into DD. From Eq.lJBJ, using T e = 
0.26 keV and assuming an integrated luminosity 
of 100 fb _1 one expects a total sample of 15000 
■0(3770) events. Assuming an integrated lumi- 
nosity of order 10 ab _1 , as suggested for future 
upgrades of B factories, more than a million of 
■0(3770) could be produced. Let us discuss a num- 
ber of potentially interesting studies. The rela- 
tive amount of charged versus neutral final states 
in the decays of -> KK, 0(3770) -> DD and 
T(45) — > BB and its energy dependence around 
the resonance is sensitive to phase space effects, 
Coulomb attraction and hadronic final state in- 
teractions. The measurement of '0(3770) in the 
radiative return and its separation into charged 
and neutral modes automatically involves the full 
energy range. 

Let us now speculate about the potential of 
such a measurement, once millions of 0(3770) are 
available 41 . Tagging the flavour of one D me- 
son through the leptonic decay of the other me- 
son, all those studies are feasible, which can also 
be investigated in the decay of neutral D mesons, 
which are tagged through the charge of the pion 
in the D* — > ixD decay. Mixing would mani- 
fest itself through like sign dilepton events from 
DD decays, and the time dependence of this phe- 
nomenon! could be studied in the present setup. 
Another option would be the study of the time de- 
pendence of doubly Cabbibo suppressed decays. 
A different subject, specific to ^(3770) decay, is 
the investigation of CP eigenstates in the decay 
of both D° and D° and the corresponding time 
evolution. Assuming CP-conservation the combi- 
nation D°(— ► tt + tt~)D°(-^ 7r + 7r~) is for example 
strictly forbidden, independent of the question of 
DD mixing. In contrast to 0(3770) production 
at a symmetric charm factory it would be possi- 
ble to study not only the rate for such a process 
but also its time dependence, which will carry 
additional information on direct and indirect CP- 
violation. Last not least one may use lepton tags 
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on one side and study the time dependence of de- 
cays into CP eigenstates at the other side, thus 
exploring combinations quite similar to those in- 
vestigated in the BB system at asymmetric B- 
meson factories. Predictions for _DZ?-mixing and 
CP violation within the Standard Model 42 lead 
to fairly small effects which would be difficult to 
observe. Nevertheless these measurements could 
lead to important limits and perhaps give access 
to physics beyond the Standard Model. 

7. CONCLUSIONS 

Measurements of the pion form factor through 
the radiative return offer the unique possibility 
for improved predictions for the muon magnetic 
moment and the electromagnetic coupling a(Mz) 
but also gives access to many new phenomena. 
The Monte Carlo event generator PHOKHARA is 
being developed further to include new features. 
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